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We study the superconducting state of S^RuCU on the bases of a phenomenological but orbital 
specific description of the electron-electron attraction and a realistic quantitative account of the 
electronic structure in the normal state. We found that a simple model which features both 'in 
plane' and 'out of plane' coupling with strengths U\\ = 40meV and U± = 48meV respectively repro- 
duced the experimentally observed power law behaviour of the low temperature specific heat C V {T), 
superfluid density n s (T) and thermal conductivity in quantitative detail. Moreover, it predicts that 
the quasi-particle spectrum on the 7 -sheet is fully gaped and the corresponding order parameter 
breaks the time reversal symmetry. We have also investigated the stability of this model to inclu- 
sion of further interaction constants in particular between orbitals contributing to the 7 sheet of 
the Fermi surface and the a and f3 sheets. We found that the predictions of the model are robust 
under such changes. Finally, we have incorporated a description of weak disorder into the model 
and explored some of its consequences. For example we demonstrated that the disorder has a more 
significant effect on the /-wave component of the order parameter than on the p-wave one. 



I. INTRODUCTION 

The symmetry of the order parameter in superconduct- 
ing Sr2Ru04 has been a subject of intense experimental 
and theoretical interest in recent years 1 ^. It is proba- 
bly the best candidate, currently, for an odd-parity, spin 
triplet, superconductor which would be a charged particle 
analogue of superfluid 3 He^ Although a number of other 
superconductors are also possible spin-triplet supercon- 
ductors (including UPt3, UGe2, ZrZn2, and Bechgaard 
salts) strontium ruthenate is probably the one which 
is best characterized experimentally. Samples can be 
grown which have exceptionally long mean free paths, 4 
and above T c the normal state is a Fermi liquid with a 
well understood Fermi surface^. 

Currently controversy exists over two key aspects of the 
Si'2Ru04 pairing state. Firstly, the gap function symme- 
try is still not known. Rice and Sigrist£ suggested sev- 
eral possible gap functions for Sr 2 Ru04 corresponding 
to analogues of superfluid phases of 3 He. Of these only 
the analogue of the Anderson-Brinkman-Morel (ABM) 
state 3 -, 



dicate that the gap function has lines of zeros on the 
Fermi surface2ii2iii. This discrepancy is not easily re- 
solved since a complete group theoretic classifications 
of all symmetry distinct pairing states in tetragonal 
crystals 1 ^* 1 ^ 1 ^* 1 ^ 1 ^ 1 ^ do not include any states which 
have both spontaneous time reversal symmetry breaking 
at T c and symmetry required line nodes on a cylindrical 
Fermi surfaced A number of 'f-wave ' gap functions have 
been propose d 19 i 20 i 21 for Sr 2 Ru04, 



d(k) ~ /(k)e z 



(2) 



d(k) ~ (k x + ik y )e z 



(1) 



is consistent with the observations of a constant a — b 
plane Knight shifti and spin susceptibility^ below T c . 
This state is also consistent with the /i-SR experiments 
which show spontaneous time reversal symmetry break- 
ing at T c m2> However this gap function has no zeros on 
the three cylindrical Fermi surface sheets^ of Sr2Ru04, 
in direct contradiction to several experiments which in- 



where /(k) is an I — 3 spherical Harmonic function. 
Such gap functions have constant a — b plane Knight shift 
and may have both time reversal symmetry breaking and 
line nodes, however in tetragonal symmetry crystals they 
are always either of mixed symmetry (requiring a double 
phase transition) or are in the same symmetry class (E u ) 
as / = 1 'p-wave' states which do not have line nodes. 
Such f-wave functions may be possible physically (de- 
pending on the details of the actual pairing interaction) , 
but the line nodes are not required by the symmetry of 
the pairing state. 

The second controversy about the Sr2Ru04 gap func- 
tion concerns the presence of three different Fermi surface 
sheets, a, p and 7. The orbital dependent superconduc- 
tivity model of Agterberg, Sigirst and Rice 2 * envisioned 
a dominant gap on one part the Fermi surface (originally 
a, /3), with the gap function on the other band only aris- 
ing from interband coupling and hence being significantly 
smaller. This theory predicted that weak impurity scat- 
tering would destroy the small gap on the inactive sheet, 
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and hence lead to a finite residual density of states at zero 
energy. However the experimental specific heat data- 
shows that Cy /T is zero at T — 0, and hence there is 
a finite order parameter on all sheets of the Fermi sur- 
face. In a recent letter, Zhitomirsky and Rica 2 ^ have 
argued that the gap function of superconducting stron- 
tium ruthenate can be described by an effective, k-space, 
interband-proximity effect. In this model they propose 
that the superconductivity is due to an attractive inter- 
action in the p-wave channel, which is acting almost en- 
tirely on one sheet of the Fermi surface, the 7 sheet. The 
other two Fermi surface sheets, a and /3 are driven to 
become superconducting because of a "proximity effect" 
or Josephson like coupling between the 7 and a,f3 bands. 
This model has a number of features which are consis- 
tent with the experimental facts, such as the presence 
of both line-nodes in the gap function and spontaneous 
time reversal symmetry breaking below T c . Furthermore, 
if the interband Josephson coupling energy is chosen to 
be sufficiently large, then the energy gap at low temper- 
atures is moderately large on all the Fermi surface sheets 
and there is no second peak below T c in the specific heat 
capacity. 

In a recent paper we have proposed a quite gen- 
eral semi-phenomenological methodology for studying 
the possible superconducting states of Sr 2 Ru04. In this 
approach one chooses, more or less systematically, or- 
bital and position dependent interaction constants to de- 
scribe the electron-electron attraction. The simplest use- 
ful model we have studied prominently featured inter- 
layer coupling 26 . This model characterizes the pairing 
interaction in terms of two nearest-neighbor negative-?/ 
Hubbard interactions, one, U\\ acts between Ru d xy in a 
single R11O2 plane, while the second, U± acts between 
Ru d xz d yz orbitals between planes. When these two pa- 
rameters are chosen so as to give a single phase transition 
temperature at the observed T c of 1.5K we find excellent 
agreement with the measured specific heat, penetration 
depth and thermal conductivity data. The gap function 
has both time reversal symmetry breaking, but also hor- 
izontal lines of nodes in the planes k z = ±7r/c on the /3 
Fermi surface sheet. The 7 sheet remains node-less, with 
a gap function of the form d(k) ~ (sinfe^ + isiak y )e z , 
corresponding to the 2-d analogue of the 3 He A-phase. 
The predicted gap function is similar to that of Zhito- 
mirsky and Rice (ZR) 25 , but differs in that it is more or 
less same size on all three Fermi surface sheets. More- 
over, while ZR rely on 'proximity coupling' to avoid the 
double phase transition we exploit the freedom provided 
by the experimental data and achieve the same end by 
fixing both Un and U± so that there is only one transition 
at the observed T c = 1.5K. 

The purpose of this paper is to clarify a number of 
unresolved questions concerning the interlayer coupling 
model. Firstly we show in Section III that the results 
of the model are quite generic, and do not depend sen- 
sitively on the choice of the specific Hubbard model pa- 
rameters which we used in Ref^&. Secondly we examine 



the effects of weak disorder on the gap function (Section 
IV) . We show that weak disorder can suppress any f-wave 
components of the gap function, while leaving the p-wave 
order parameter relatively unchanged. Finally in Section 
V we study the generalisation of our model by allow- 
ing for a "bond proximity" interactions. Such symmetry 
mixing interactions have been proposed by Zhitomirsky 
and Rice 2 ^ as a mechanism leading to the single super- 
conducting transition temperature. It turns out that the 
mechanism operates in the orbital picture as well and we 
obtained the single superconducting transition tempera- 
ture, but different slope and jump of the specific heat. 

II. GAP SYMMETRY AND PAIRING BASIS 
FUNCTIONS 

Let us begin by reviewing briefly the symmetry prin- 
ciples which are used to classify different pairing symme- 
try states in odd-parity superconductors. We shall use 
these principles to contrast the different pairing states 
that have been proposed for strontium ruthenate. 

On very general ground we expect that the phase tran- 
sition into the superconducting state is of second order, 
and so there exists an order parameter, or set of or- 
der parameters, r]i(r),i — l,...n. For superconductors 
these order parameters are complex, transforming under 
the U(l) gauge symmetry as r\i — > e r\i. Therefore the 
Ginzburg-Landau Free energy can always be expanded 
as 

F s = F n + / d 3 r c^*(r)d fe 77,(r) 

J \2m i: jki 

+ f) ijkl r ] *(T)r l *(r)r, k (r)r lk (r) + ...) (3) 

where summation convention is implied for the indices i,j 
etc, and as usual di = V ' i — 2eiAi/Ti, with A the magnetic 
vector potential. 

If the normal state above T c possesses a symmetry 
group Q, then the order parameters rji can be grouped 
into terms corresponding to the different irreducible rep- 
resentations r of Q, transforming under symmetry oper- 
ations as 

'/! •"!.;;;/"/; (4) 

where g £ G, and the matrices R^j(g) constitute the rep- 
resentation r of the group Q. 

The general theory of group representations implies 
that we can choose a basis in which the matrix ol%j is block 
diagonal, with each block corresponding to an irreducible 
representation, T. In this basis the full Ginzburg-Landau 
Free energy is of the form 

F s = F n+ [ d*r f^^^f (r)<W(r) 
J \r,r> zm ijki 
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TABLE I: Irreducible representations of even and odd parity 
in a tetragonal crystal. The symbols X, Y Z represent any 
functions transforming as x, y and z under crystal point group 
operations, while / represents any function which is invariant 
under all point group symmetries. 



Rep. 


symmetry 


Rep. 


symmetry 


A 29 
Bi 9 


I 

XY(X 2 -Y 2 ) 
X 2 - Y 2 


A iu 
A 2u 
Biu 


XYZ{X 2 - Y 2 ) 
Z 

XYZ 


B 2s 
E y 


XY 

{XZ, YZ} 


B 2u 

E u 


Z(X 2 - Y 2 ) 
{X,Y} 



TABLE III: Basis functions ■yf (k) for the odd parity irre- 
ducible representations of body-centred tetragonal crystals. 

Rep. in-plane inter-plane 
Ai„ 

A 2u - cos if cos ^- sin i^- 

B lu - sin if sin ^ sin ^ 

B 2u 

E u sin k x sin if cos cos if^ 

sin k y cos if sin ^ cos if^ 



TABLE II: Products of the irreducible representations of D&h 
point group symmetry 



® 


Ai 


A 2 


Bi 


B 2 




E 


Ai 


Ai 


A 2 


Bi 


B 2 




E 


A 2 


A 2 


Ai 


B 2 


Bi 




E 


Bi 


Bi 


B 2 


Ai 


A 2 




E 


B 2 


B 2 


Bi 


A 2 


Ai 




E 


E 


E 


E 


E 


E 


Ax € 


B A 2 © Bi © B 2 



+ ^<.r ? r(r)r ?3 r (r) 
r 

rrT'T'" / 

(5) 

The quadratic term involves only a single represen- 
tation, r. At T c , in general, only a single irreducible 
representation will have a zero eigenvalue of the block 
diagonal matrix oq-. Therefore only the components of 
the order parameter Tjf corresponding to that eigenvector 
will become non-zero just below T c . 

Now let us apply these very general principles to the 
specific case of spin-triplet pairing in Sr2Ru04. This is a 
body-centred tetragonal crystal with inversion symmetry. 
The relevant crystal group is D±h, and Table |U shows its 
irreducible representations. For each representation its 
symmetry is denoted by a typical function, where the 
symbols X, Y, Z represent any functions which trans- 
form as x, y and z under the point group operations, 
and I means any function which is invariant under all 
point group operations. The representations A\ g . . . E g 
have even parity, while Ai u . . . E u have odd parity. Table 
ITT1 shows the multiplication table for the irreducible rep- 
resentations, i.e. how direct products of representation 
matrices r (g) F' decompose into a sums of block diagonal 
matrices Ti © T 2 © . . .. 

An immediate consequence of the multiplication table 
HTlis that in tetragonal crystals the order parameter is ei- 
ther of a single representation T only, or there are two or 
more distinct thermodynamic phase transitions. This is 
because to quadratic (or higher) order in the Ginzburg- 
Landau free energy there are no symmetry allowed cou- 
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FIG. 1: Body-centred tetragonal lattice, showing the nearest 
neighbour pairs in-plane, and between planes. 

pling terms of the form 

^ rr ^'*(r)^(r) % r (r) ?7 r(r) 

in Eq. [SJ The prooiii is simply that r' <g) T ® V <g) T never 
contains the identity representation Ai g , and hence such 
terms are not allowed as quartic invariants of the Free 
energy (or at higher order). In the absence of such terms 
the free energy functional is always of at least quadratic 
order in the subdominant order parameter rjf (r), and 
hence these subdominant components can only become 
non-zero in a separate phase transition below T c . 

Using these irreducible representations we can expand 
the BCS gap function in terms of functions of each sep- 
arate symmetry class. For odd parity pairing states we 
can represent the BCS gap function by a vector d(k) or 
a symmetric 2x2 complex matrix 

/ A TT (k) A u (k) \ / idyQt) ~ d x (k) d,(k) \ 

^ A U (k) A U (k) d,(k) d x (k) +l dy(k) J 

(6) 

where A n (k) = A iT (k) and A (T(T /(k) = -A CT(r '(-k). For 
each irreducible representation we can choose a complete 
set of orthonormal basis functions in the Brillouin zone, 
7 2 r (k). Expanding the gap function in terms of these 
functions we have 

A^(k)=X>Lr,<yf(k). (7) 

i 

The expansion coefficients essentially provide the set of 
order parameters in Eq[SJ The basis functions must be 
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periodic in reciprocal space, 7 4 r (k) = jf (k+G), or equiv- 
alently, they must obey periodic boundary conditions in 
the 1st Brillouin zone. They can be chosen, most nat- 
urally, in terms of their real-space Fourier transforms, 
which correspond to lattice sums of the real-space Bra- 
vais lattice. For a body-centred tetragonal crystal, such 
as Sr2Ru04 shown in Fig. ^ the leading basis functions 
correspond to the four nearest-neighbour in-plane lattice 
vectors, R = ±ae x and R = iae v , giving two odd par- 
ity basis functions: sin/c x a and sinfcj,a. The eight body- 
centred lattice vectors R = ±%e x ± |e y ± |e 2 lead to the 
four odd-parity basis functions shown in the last column 
of Table ITTT1 (where for simplicity we have chosen units of 
length such that a = 1 ) . In the models which we investi- 
gate in the remainder of this paper, we shall assume that 
these basis functions, Table ITTT1 are sufficient to describe 
the gap function. Physically this corresponds to the as- 
sumption that the paring interaction V aa -' (r, r') is short 
ranged in real-space. 

Considering Table [I] we can see that in S^RuO^ "p- 
wave" pairing states can correspond to either the A2„, 
(or p z ) representation or the doubly degenerate E u repre- 
sentation (p x , p y ). The only symmetry distinct "f-wave" 
pairing states are the Bi„ and B2 U representations, cor- 
responding to f xyz and f( x 2- y 2) z type symmetries. Nei- 
ther of these states can be used in the case of a two- 
dimensional single-plane model of S^RuO^, since they 
both become zero in the plane k z = 0. It is also interest- 
ing to note that in Table ITTT1 there are no basis functions 
of Ai u or Bi u symmetry. Pairing in these channels would 
require long range interactions extending to at least the 
inter-plane second nearest neighbors. 

In the light of these symmetry principles let us com- 
ment on a number of the possible gap functions which 
have been proposed for Sr 2 Ru04. Among the five states 
described by Rice and Sigrislpi the only one consistent 
with the Knight shift experiments isi£ 

d(k) = (sinfca; + i sin k y )e z (8) 

belonging to the E u representation of Table ITTT1 It breaks 
time reversal symmetry, consistent with the p-SK exper- 
iments of Luke et al*^, and leads to a spin susceptibil- 
ity which is constant below T c for fields in the a — b 
plane, consistent with Knight shift 7 and neutron scat- 
tering experiments 8 . However it has no gap nodes on 
a Fermi surface of cylindrical topology, such as the a, 
(3 and 7 sheets of Sr2Ru04, and therefore is inconsistent 
with the heat capacity^ penetration depthlS and thermal 
conductivity experiments^. 

On the other hand the f-wave gap function proposed 
by Won and Maki^fi 

d(k) ~ k z (k x ±ik y ) 2 e z (9) 

has both line nodes and broken time reversal symmetry 
below T c . However from the symmetry analysis above, 
it is clear that this does not correspond to a single irre- 
ducible representation of the symmetry group. It is a sum 



of the function k z (k x — fc 2 ,), belonging to B2 U and k x k y k z 
belonging to Bi„. Although they would be degenerate 
in a system with cylindrical symmetry, in a tetragonal 
crystal they will be non-degenerate and hence have dif- 
ferent T c s. The Bi„, B2 U states individually possess time 
reversal symmetry. Therefore with this order parameter 
we would expect to find a specific heat anomaly with two 
transitions, and time reversal symmetry breaking would 
only occur at temperatures below the lower transition. 

The f-wave order parameter proposed by Graf and 
BalatskyiS, 

d(k) ~ k x k y (k x + ik y )e z (10) 

is in the same symmetry class as E„, since B2 <S> E = E 
in Table HTl Therefore in the sense of pure symmetry ar- 
guments the gap nodes in planes k x = and k y = are 
"accidental" . Such a gap function is certainly valid, but 
the nodes are present for reasons connected with the spe- 
cific microscopic pairing interaction employed, and not 
required by symmetry alone. This comment also applies 
to the Bi eg) E f-wave state 

d(k) ~ {k 2 x - k 2 y )(k x + ik y )e z (11) 
discussed by Dahm, Won and Makij2I and Eremin et 

The full group theoretic classification in tetragonal 
crystalsi2*i^*i^iSiiLi 8 i and the above analysis does not 
show a single pairing state with both symmetry required 
lines of nodes and spontaneously broken time reversal 
symmetry below T c . Therefore, if we accept both the /i- 
SR and low temperature thermodynamic and transport 
measurements, then we must consider states which have 
lines of nodes for specific microscopic reasons, rather then 
for pure symmetry reasons. 

In the remainder of this paper we shall focus on the 
specific model which we proposed in a previous paper—, 
in which the lines of nodes appear in the plane k z — 
±7r/c, derived from the pair of inter-plane basis functions 
of E u : 

• ky k^C . Jv^C 

sin — cos — cos , cos — sin — cos 

2 2 2 2 2 2 

from Tabl dIIl| as originally suggested by Hasegawa et 
alSL. 



III. INTERLAYER COUPLING HAMILTONIAN 

Since the underlying microscopic mechanism for su- 
perconductivity in Sr2Ru04 is not known we choose to 
adopt a phenomenological approach to the pairing mech- 
anism. We first make an accurate tight binding fit to the 
experimentally determined Fermi surfac e) 5 ! 28 and then in- 
troduce model attractive interactions between the differ- 
ent orbitals centered on different sites. We can inves- 
tigate different 'scenarios' depending upon which model 
interactions are assumed to dominate. Frequently, when 
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symmetry to have the following form 




FIG. 2: The Fermi surface of Sr2Ru04 in the plane k z = 0, 
obtained by fitting the de Hass data of Bergman et al. 5 . Note 
that the alpha Fermi surface sheet has only two-fold symme- 
try, because of the shape of the Brillouin zone boundary. 



these pairing interaction parameters are chosen to repro- 
duce the experimental T c , there is no freedom to adjust 
the parameters further. Once the parameters have been 
selected, then a number of different experimental quan- 
tities can be calculated independently and compared to 
experiment. The goal is to find one specific paring sce- 
nario which agrees with all of the experimental obser- 
vations. If this can be achieved then one has found an 
effective Hamiltonian for the pairing, which can be inter- 
preted physically. This effective pairing Hamiltonian can 
then be used to guide the search for the true microscopic 
Hamiltonian. This methodology has proved very useful 
in cuprate superconductivity 29 and here we shall deploy 
it to study Sr2RuC>4. 

The effective pairing Hamiltonian we consider is a sim- 
ple multi-band attractive U Hubbard model: 



H = E (( £ m - l^)SijS n 
ijmm' ,<r 

-\ E 

ijmm' era' 



tmm'(ij)) Cimo-Cjm'a 



(12) 



where m and ml refer to the three Ruthenium t2 g orbitals 
a = xz, b = yz and c = xy and i and j label the sites of 
a body centered tetragonal lattice. 

The hopping integrals t mm >(ij) and site energies e m 
were fitted to reproduce the experimentally determined 
Fermi Surfac o 5 i 28 . The nearest neighbour in-plane hop- 
ping integrals along R = e x , where the ab plane lattice 
constant is taken to be 1, are constrained by the orbital 




(13) 



(and similarly for R = e y taking into account sign 
changes due to orbital symmetries). The next nearest 
neighbour in-plane hopping integrals along e x + e y were 
assumed to be of the form 



t ab 

= I U 
t' 



(14) 



The parameter t' controls the shape of the 7-band Fermi 
surface, while the parameter t ab determines the hy- 
bridization between the a and b orbitals and hence the 
shape of the a and (3 Fermi surfaces. The c — axis mag- 
netic field de Hass van Alphen data™ gives the areas and 
cyclotron masses of the three Fermi surface sheets, and 
these six numbers can be fit exactly with t — 0.08162eV, 
t' = -0.45*, t ax = lMt, t bx = 0.06t ax , t ab = 0.08t ax , 
and the on-site energies were e c = — 1.615t and e a = 
e b = -1.062t QX . 

To obtain a three dimensional Fermi surface we as- 
sumed that the dominant inter-plane hopping is along 
the body-centre vector R = \{e x + e y +ce z ) and has the 
form 



£_L tfiyb thyb 
thyb tj- thyb 
thyb thyb 



(15) 



and similarly for R = ^(±e x ±e y ±ce z ) with appropriate 
sign changes. The parameter thyb is the only term in 
the Hamiltonian which mixes the c orbitals with a and 
b. With only these two parameters it is not possible 
to fit exactly the full three dimensional Fermi surface 
cylinder corrugations determined by Bergemann et alA, 
but the parameters thyb = 0.12i o &, t± = — 0.03t a {, give 
a reasonable agreement for the dominant experimental 
corrugations. Fig.|21shows the fitted Fermi surface in the 
plane k z = in the extended zone-scheme. Note that the 
a sheet has only two-fold symmetry, due to its position 
centred on the Brillouin zone boundary at X. 

The set of interaction constants U™ m ,{ij) describe at- 
traction between electrons on nearest neighbour sites 
with spins a and a' and in orbitals m and w! . Thus our 
actual calculations consists of solving, self-consistently, 
the following Bogoliubov-de Gennes equation: 



E 

jm' a' 



rpu TT 



m.m 

Hrn.rn, 



(ij) 



jm' a' 
jm' a 1 



o. 



(16) 

where H mm i(ij) is the normal spin independent part of 
the Hamiltonian, and the &.™ m ,(ij) is self consistently 
given in terms of the pairing amplitude, or order param- 



(17) 
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defined by the usual relation 



A mm 



(ij) 



EU I/* 
^imcr^jra 



,A±-V{e v )). 



(18) 



where v enumerates the solutions of Eq. 1161 

We solved the above system of Bogoliubov de Gennes 
equations including all three bands and the three di- 
mensional tight-binding Fermi surface. We considered 
a large number of different scenarios for the interaction 
constants. First we assumed that the pairing interaction 
Umm' {ij) for nearest neighbours in plane is only acting 
for the c (d xy ) Ru orbitals. In this case both a d-wave 
(d x 2_ y 2) pairing state and p-wave ((k x + ik y )e z ) states 
are possible. The d-wave state has line nodes, but would 
not be consistent with the experiments showing constant 
Knight shift and time reversal symmetry breaking be- 
low T c . Therefore we discard such solutions here, and 
only concentrate on the odd-parity spin triplet solutions. 
The motivation is not to explain the microscopic pair- 
ing mechanism, but to model pairing state produced by 
various types of effective attractive interactions. These 
attractive interactions may arise from, for instance, fer- 
romagnetic spin fluctuations&SSiSiSl, which can favour 
spin triplet pairing compared to the d-wave solutions. 
However, their origin may be more complicated, for ex- 
ample a combined electron-phonon and spin fluctuation 
mechanism. 

With only the nearest neighbor in-plane interactions 
the set of possible odd-parity, spin triplet, solutions that 
we found never includes any possible state with nodes 
of the gap. Therefore we extended the model to include 
inter-plane interactions. Using two interactions, a near- 
est neighbor in-plane interaction, (i — j in Fig. 0, and a 
nearest neighbor inter-plane interaction, (i — ! in Fig. ^) 
which fulfill the tetragonal symmetry, we have the two 
types of basis functions for the gap equation given in Ta- 
ble IIIII Then we have the possibility of horizontal line 
nodes in the gap arising from the zeros of cos (k z c/2) at 
k z = ir/c on a cylindrical Fermi surfaced. 

Because the pairing interactions U™ m ,{ij) were as- 
sumed to act only for nearest neighbor sites in or out of 
plane, the pairing potential A™ m ,(ij) is also restricted 
to nearest neighbors. We further focus on only odd 
parity (spin triplet) pairing states for which the vec- 
tor d ~ (0,0, d'), i.e. A^ m ,(ij) = A£ w ,(y), and 
Aj] m ,(y) = A~ m ,(ij) = 0. Therefore in general we 
have the following non-zero order parameters (i) for in 



plane bonds: A^ m , (e x ), A^ lm ,(e y ), and (ii) for inter- 
plane bonds: A^ m ,(Ry) for Ry = (±a/2, ±a/2, ±c/2). 

Taking the lattice Fourier transform of Eq.^lthe cor- 
responding pairing potentials in k-space have the general 
form (suppressing the spin indices for clarity): 
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Al^,sinfc, + Al p «,sinfc y 



+A 
+A 



mm' x 1 mm' 



mm' sm Y cos y cos Y 



k z c 
2 
k 7 c 



vara' 



k z c 

2 
k. 



i / ky kg C / \ 

- A mm' Sm Y Sm 2 Sm lT"' ( ' 

Note that beyond the usual p-wave symmetry of the 
sin k x and sin k y type for the c orbitals, we include all 
three additional p-wave symmetries of the sin/c/2 type 
which arc induced by the effective attractive interactions 
between carriers on the neighboring out-of-plane Ru or- 
bitals. These interactions are also responsible for the 
f-wave symmetry order parameters, A^f m ,, transform- 
ing as B lu in Table 1. This latter is symmetry dis- 
tinct from all p-wave order parameters in a tetragonal 
crystal, unlike the other f-wave states discussed in the 
introductio n 19 ! 20 ! 21 ! 22 . The p z order parameters A m ^, 
are of A 2u symmetry. In contrast the pairs A x ^, , A^, 

are of the same E u 'p-wave' symmetry as Ajj^, , Aj^, . 
In general, the order parameters in each distinct irre- 
ducible representations have different transition temper- 
atures, as expected from Eq. [S] 

In a recent paper— we chose a particularly simple set 
of attractive pairing interactions U™ m ,(if). For in-plane 
nearest neighbours we assumed that the pairing interac- 
tion is only acting for the c (d xy ) Ru orbitals only 



Un 






U n 



where U\\ = 0.494i, 



(20) 



On the other hand, given that the ruthenium a and b 
orbitals (d xz ,d yz ) are oriented perpendicularly to the 
planes we choose to introduce the inter-plane interaction 
only for these orbitals, 



U± U±. 
U± U± 




where U±_ = 0.590t. (21) 



„„„, smy COSy cos — 



Therefore we have, as a minimal set, only two coupling 
constants U\\ and U± describing these two physically dif- 
ferent interactions. 

As discussed earlier our strategy is to adjust these phe- 
nomenological parameters in order to obtain one tran- 
sition at the experimentally determined T c . Thus, be- 
yond fitting T c , there are no further adjustable parame- 
ters, and one can compare directly the calculated physical 
properties of the superconducting states to those exper- 
imentally observed. Consequently, if one obtains a good 
overall agreement one can say that one has empirically 
determined the form of the pairing interaction in a phys- 
ically transparent manner. Evidently such conclusion is 
the principle aim of the calculations. 

As we have shown in Ref.26, this two parameter 
scenario gives an excellent agreement with the experi- 
mental specific hea1i£, superfluid density^ and thermal 
conductivity— We chose the constants U\\ and U±, so 
that there is a single phase transition at T c = 1.5K, corre- 
sponding to the values given in Eas.l20landl21l Below T c 
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x 10"' 
0.5, — 



x 10"' 




0.33 



FIG. 3: (a) Temperature dependence of order parameters 
|Aaa x |, |Aj!c x | and |A^. Pa= | for a number of u values (u' = 
u). (b) Order parameters lAU^I and lAatf" 1 !, |A^ c Px | at zero 
temperature versus the interaction parameter u(= it'). 



the order parameters have the symmetries A 



\\p« 



A^ b Py = iA^l x as expected for an E u pairing symmetry«2i 
(k x + ik y )e z corresponding to the same time reversal bro- 
ken pairing state as 3 He — A. We also found that a 
much lower temperatures, additional transitions occurred 
where the /—wave and p z order parameters become non- 
zero. The gap function has line nodes on the Fermi 
surface, in agreement with experiment, only when the 
f-wave component is zero. Arguing that the f-wave com- 
ponent would be suppressed by impurities, we showed 
that with the f-wave component removed, one obtains 
excellent agreement between the calculated and experi- 
mental specific heat, penetration depth and thermal con- 
ductivity. We show, in Sec. IV below, that this removal 
of the f-wave component is justified by the presence of 
weak disorder. 

It is important to ask how these results depend on the 
details of the assumptions made in the model. In order 
to test the stability of our results to variations in the 
model we therefore introduced some additional subdom- 
inant interaction parameters. For our initial exploration 
of the issues involved we have generalized Eqs. I|2U[) and 
(|2*T|) as follows: 




(22) 



U± U± v! 
U± U± u' 
v! v! u' 



Reassuringly, with these modified parameters we ob- 
tained a temperature dependence of the gap parameters 
which are qualitatively similar to those for the original 
parameters. It is interesting to note that for fixed values 
of U± and Un the changes of u and v! hardly change the 
superconducting transition temperature. We have sys- 
tematically studied the effect of additional interactions, 
especially so on the line u = u', and found small differ- 
ences compared to the u — solution even for u as large 
as 0.28i . The differences are mainly connected with the 
appearance of out of plane components of generated 
by the new interactions as is evident from Fig.©. For 
larger values of u the difference becomes more signifi- 
cant (Fig. |2K). Note, however, that only low temper- 
ature dependence of the pairing amplitudes is affected. 
In Fig. [3p we show the variation of a few characteris- 
tic |A mm | against u at zero temperature. Clearly, for 
u > 0.3t there is a qualitative change of our solution 
leading to dominant out of plane pairing components in 
all orbitals. Large u also affects the critical temperature 
T c . Interestingly, for finite u we also observe increasing 
values of in-plane pairing amplitudes in the a and b chan- 
nels: A^ p ^, and A^ p ^, for to, to' = a, b. Reassuringly, 
the corresponding specific heat (Fig. is essentially un- 
changed and remains in equally good agreement with the 
experiments. Therefore we conclude that the solution we 
have found is not very specific to the precise details of 
the model parameters which we assumed, but is a generic 
solution valid for at least some range of the possible in- 
teraction parameters of the form depicted in Eq. 1221 

The quasiparticle energy gap structure which we ob- 
tained is shown in Fig. [S] The gap is finite everywhere 
on the 7 sheet, Fig.|^d), although it is very anisotropic, 
and becomes small when the Fermi surface approaches 
near to the van Hove points at (jr, 0) and (0, 7r). In con- 
trast, the a and (3 Fermi surface sheets have gap zeros in 
the vicinity of the lines k z — ±7r/c. In the case of f3 the 



^70 
O60 

E 

^40 



-20 



O 





u/t=0.0 

u/t=0.28---- 

. u/t-0.30 

v. . . 







0.0 



0.4 



0.8 

T/Tc 



1.2 



1.6 



FIG. 4: Calculated specific heat for a three parameters u(= 
u) (u/t — 0.0, 0.28 and 0.30 corresponding to full dashed 
and dotted lines, respectively) compared to the experimental 
data (points) of NishiZaki el al..[8] 
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a -band 
(a) 



a -band 
(b) 



/ r 



P -band 
(C) 



y-band 
(d) 




x 10 
0.6, 



kz= rfc+0.150 .. 
kz= rfc+0.085 - 
n kz= it/c - 




1 "° ^ 2 '° 



FIG. 5: Lowest energy eigenvalues, E v (k) on the Fermi 
surface; a sheet in the plane k z = tt/c + 0.085 (a) and 
fc 2 = tt/c + 0.085 (b), /3 (c) and 7 (d) sheets in the plane 
k z = 0. 



line fc z = —tt/c. This interesting nodal structure arises 
from the fact that the a Fermi surface cylinder is cen- 
tered at X in the Brillouin zone not at T (Fig. and 
therefore it has two-fold symmetry not four fold like j3 
and 7. Notice also that the excitation gap on the a sheet 
is non-zero even when A aQ = A a fc = A» = 0, because it 
is hybridized to the c orbital and A cc 7^ 0. 





Ey! .-- ^^^^^^W 
kz " ^^^^^^^M 



( c ),# 



FIG. 6: Minimum energy quasiparticle eigenvalues on the 7 
Fermi surface sheet, E(k.w), plotted in cylindrical polar co- 
ordinates as functions of k z and a — b plane polar angle, 6. 
Parameter values are u/t = 0.28 (a), 0.30 (b), 0.32 (c), re- 
spectively. One can see that for u < 0.3 the 7 sheet gap is 
nodeless, while for u > 0.3 line nodes appear at k z — ±7r/c. 



Note that this nodal structure of the gap is unchanged 
by the presence of the small subdominant interaction pa- 
rameter u, in Eq. l22|) . However, upon increasing the 
value of the u parameter eventually the results change 
qualitatively, leading to appearance of additional line 
nodes in 7 (Fig. iJBJ) for u = 0.322;. In this case the 
7 band gap also develops a line node, similar to the be- 
havior of the band. 



IV. EFFECTS OF DISORDER 



gap is zero to numerical accuracy on these nodal lines. 
While in the case of a the gap is very small on these lines, 
but not exactly zero. In fact there are eight point nodes 
on the a sheet, as can be seen in Fig.|^a,b). Two point 
nodes lie just above the k z —tt/c line at k z w tt/c + 0.085 
at two different angles. Another pair lie just below, at 
k z w tt/c — 0.085 at an angle rotated by <f> = tt/2. The 
remaining four are located in similar positions near the 



As we noted it earlier, to obtain agreement with ex- 
periment we had to eliminate the /-wave component 
A^-^,(T) and we suggested that this can be done by an 
appeal to the effects of a small amount of disorder. We 
shall now substantiate this contention by explicit calcu- 
lations. 

In case of non-magnetic disorder our Hamiltonian can 
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FIG. 7: Order parameters |A^ a P:c |, \A l l^\, |A^/| for the clean FIG. 8: Zero temperature order parameters |A^f*|, |Al?*|, 
and disordered systems (t" 1 = 0.005) as functions of temper- |A„/| as a function of impurity scattering rate r _1 . 
ature. 



be written 

H. ^ ((^m ~t~ L^&ij&mm' ^mm'(u)) ^ima^"j m ' a 

ijmm' ,cr 



\ E U m m >(i3)f> 



%ma ' L jm' a' 



(23) 



ijmm' aa' 



where q is a random site energy. For a given configu- 
ration of €i one can, in principle, perform calculations 
(Eq. I17I18|I and then average over many configurations. 
More readily, for highly disordered systems it is possible 
to apply mean field theory of disorder by making use of 
the Coherent Potential Approximation CPA 31 i 32 i 33 i 34 . 

Here however, as superconducting S^RuCU samples 
were found to be relatively clean, we can limit our anal- 
ysis to weak disorder and non-resonant impurity scatter- 
ing. Then knowing the scattering rate r _1 we can apply 
the Born approximation 3 ^ in calculating the self-energy 
of the disorder averaged Green function. Following this 
Abrikosov-Gorkov approach, we assume that impurity 
scattering will create a finite imaginary self energy of the 
order 



£(k<j) = it 1 sgn(a 



(24) 



Thus our equation of self-consistency in the configura- 
tionally averaged pair potential can be written in terms 
of Matsubara frequencies uj n — (n/ 0)(2n + 1) as follows 



x -E- 

^ i 



n— — oo 



(25) 



^%m<y u jm a' 



^iK + T-Vkl)-^' 

where S denotes a positive infinitesimal. Exchanging the 
summations over v and n indices Eq. 1251 can be writ- 
ten, after decoupling the summation over negative and 
positive Matsubara frequencies, as 



1^ 

'0 



(26) 



E 



E" 



^(a,„+r-i) 2 + W 

n— 



Conveniently the sum on the right hand site of Eq. 
can be evaluated 3 ^ and it leads to the final formula: 



%m<7 jm' a' 



(27) 



1/1 \E V 
x — Im* ' 



2?r V2 2ttt 2vr 



Note that in the limit of a clean system r 1 — » 



J_ M f^ + J_ + hEl) 2 / ( E' i i (2s, 

Z7T \ 2 Z7TT Z7T 



and Eq. |2HI coincides with that of the clean system (Eqs. 

HUB- 

Thus, in a weak disorder limit we have again solved 
the Bogoliubov-de Gennes equations including a small 
t _1 . The results of our calculations are shown in Fig. 
where we have plotted the order parameters |A„*' a! |, 

|Alc x |, |A^/| versus temperature. Evidently in the dis- 
ordered case the small /-wave amplitude (|A^/|) is re- 
duced to zero much more rapidly than the larger p-wave 
ones ( |A^|, |Al? x |). Furthermore, in Fig. H we show 
disorder dependence of these three paring amplitudes at 
zero temperature. Clearly each order parameter is re- 
duced to zero in a typical Abrikosov-Gorkov like manner, 
becoming zero approximately when the pair-breaking pa- 
rameter r|A dean (0)|/7r 1. 

From Figs. Q and © it is clear that, for moderate 
scattering rates, there is a region where the f-wave gap 
components are reduced to zero but the larger p-wave 
components are more or less unaffected. Thus we con- 
clude that the simultaneous neglect of A^^, and A^{ n , 
and t -1 is justified^. It would be an interesting experi- 
mental confirmation of this model, if ultra-clean samples 
were found to have a second phase transition at a much 
lower temperature than T c ~ 1.5K. 
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BOND-PROXIMITY EFFECTS 



x 10"' 
0.4,_ 



To get a single superconducting transition tempera- 
tures within our interlayer coupling model we are forced 
to fine tune two interaction parameters. However, it has 
been proposed 25 that a single transition can be obtained 
in an multiband model by allowing for a symmetry mix- 
ing interaction of the type J7(k, q) = g' f(k)g(q), where 
/(k) and g(q) are order parameter symmetry functions 
for respective bands. 

It is the aim of the present section to check to which 
extend similar approach may be used in our bond model. 
We start with short discussion of the source and magni- 
tude of symmetry mixing interaction. The description we 
have used is a real space, two point near neighbour inter- 
action such as naturally arises in any multi-orbital, ex- 
tended, negative U Hubbard model, Eq. To be quite 
clear about this matter we recall that a generic pair-wise 
interaction like C/(r, r'), when expressed in the language 
of a tight-binding model Hamiltonian will, in general, 
give rise to four point interaction parameters Uij^i- The 
original Hubbard Hamiltonian makes use of the one point 
parameters — Una whilst the extended Hubbard 

model is based on two point parameters U\ ? = Uij^j. 
Evidently our 'bond' model is a negative U-version of 
the latter—. The symmetry mixing interactions 2 ^ arise 

(3) 

from 3-site interactions U i ; ■ \. The physics of this is of- 
ten referred to as assisted hopping^. If one assumes, 

as is normally the case in an isotropic substance, that 
|[/(i)| > |[/(2)| > |{7(3) | > ([7(4) | then the < bonds > repre _ 

sent stronger coupling than assisted hopping and should 
be the preferred coupling mechanism. However, for the 
tetragonal arrangement of Ru atoms in Sr2Ru04 this is 
no more than a suggestion at present. 

(3) 

In the presence of a three point interaction U\^ = 
Uij,a — Ui, for all nearest neighbours ijl such that i and 
j are in one Ru plane while I is on a neigbouring one 
(Fig. nj, the gap equation (Eqs. llblTHI) can be rewritten 
in k-space as^S 

q 

+ E Cm', o' (q, k - <i)xZ' (q) • (29) 

q,oo' 

where, as before, 

X ™ m , (k) = < m X,W (i - V{E V ) . (30) 

In a body centered tetragonal crystal (Fig.QJ the vari- 
ous matrix elements of the general four point interaction 
U m m',oo' responsible for p-wave paring can be written 
(suppressing spin indices for clarity): 

E/ cc (k,q) = 2f/||V(k)V(q) (31) 
U mm > (k, q) = 8U±V(k)V(q) for m, m' = a, b 




T/t xlfj 3 

FIG. 9: Order parameters \A^*\, |A^|, for U± = 0.400t, 
Uu = 0.494£ and the proximity coupling as functions of tem- 
perature. The three curves '1', '2' and '3' correspond to the 
values 0.0, 0.005, 0.010 of three point coupling constant. 

U mm \cc{K q) = 8UiV(k)V(q) for m', m = a,b 
J7 cc , mm '(k,q) = 8L r /I/(k)V r (q) for m',m = a,b 

where V(k) and V(k) are respectively: 

V(k) = (sin k x + sin k y ) (32) 

f>/l \ ( • k% , . ky ky \ k Z C 

V (k) = sin — cos — + sin — cos — cos . 

v ' \ 2 2 2 2 J 2 

Note that the three point interaction leads to an extra 
interlayer coupling proportional to Uj. Interestingly, the 
general form of the order parameter is the same as previ- 
ously derived (Eqs. (|16I19[0 despite the additional three 
point coupling Eq. I.'i2[| in the self-consistency relation 
Eqs. UU). 

It has to be noted that the presence of interaction 
Ui strongly changes T c . To get its correct value (1.5K) 
for the present model we have taken U± = 32meV and 
U± = 40meV and repeated our calculations for various 
Ui values. The results are shown in Figs. 10 and (|10|) . 

Fig. © shows the results for the amplitudes Acc x (T) 
and A^f x (T) including the three-point interaction. Note 
that for Ui — (curves labeled by (1) in the figure) the 




T/Tc 

FIG. 10: Specific heat and the effect three point interac- 
tion, compared to the experimental data of NishiZaki et al. 
[8]. The arrow indicates increasing values of the proximity 
coupling ('1-2-3' as in Fig. 0. 
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x 10" 2 
0.4. 



0.3 




T/t x10" 3 

FIG. 11: Order parameters |A^|, |A|?*|, for U± = 0.4004, 
Uu = 0.494t and the proximity coupling as functions of tem- 
perature in disordered system r _1 = 0.0054 The three curves 
T', '2' and '3' correspond to the values Ui/t: 0.0, 0.005, 0.010 
of three point coupling constant. 



temperature where A C c (T) becomes non-zero is much 
higher than that where A^* (T) becomes non-zero. It is 
evident from the figure that for Ui ^ the parameters 
A^ m , (T) for a, b and c orbitals vanish at the same tem- 
perature and two transitions merge into one. Thus, the 
proximity coupling mechanism identified by ZR 2 — in their 
band description of the electron-electron interaction, also 
works in our bond model. 

Intriguingly, although the proximity coupling works in 
principle, the above mechanism does not seem to be help- 
ful in the context of building phenomenological interac- 
tions suitable to describe experimental data. To illustrate 
this point we reproduce, in Fig. 1)11) [I. the specific heat 

corresponding to the set of Ale" (T) and A^f 1 (T) shown 
in Fig. 10. Clearly, as Uj increases the second transition 
at low temperature becomes enlarged and merges with 
the first transition at higher temperature. However, the 
small values of Ui shown in Figs. l(51 ll0| l are not sufficient 
to get the specific heat jump at T c right. Therefore Ui 
must be large for the ZR scenario to fit the experiments. 
By contrast, as we have demonstrated earlier, if U% = 
and the sizes of U± and Uu are adjusted so that only one 
transition occurs both the low temperature slope and the 
jump at T c agrees with experiments. Thus although we 
have not investigated models featuring a 'proximity ef- 
fect' induced by type of interactions systematically 
we conclude that such interactions are not needed to fit 
the available data. 

Finally, it is also interesting to see what are the effects 
of disorder on the 'orbital proximity effect'; the results 
are shown in Fig. ^2 We see that, disorder can eliminate 
the gap on a, [3 sheets, while leaving it almost unchanged 
on 7. This feature of the proximity effect scenario opens 
it up for experimental verification by measurement on 
samples with increasing disorder. Evidently the effect of 
disorder should be that the low temperature power laws 
disappear due to the destruction of superconductivity on 
the a, (3 sheets. 



VI. Conclusions 



We have introduced a methodology for building semi- 
phenomenological, attractive electron-electron interac- 
tions bond by bond for calculating superconducting prop- 
erties under circumstances when the physical mechanism 
of pairing is not known. We deployed it to study p-wave 
pairing in S^RuO.^. A bond was described by an inter- 
action constant U^ m ,(ij) which depends on the sites i 
and j, the orbitals m and m' , and their spin orientation 
a and a'. We have solved the appropriate Bogoliubov 
de Gennes equations for a number of scenarios defined 
by a small set of interaction constants. We have found 
that the one for which Uj}(ij) = U± for i and j being 
nearest neighbour Ruthenium atoms in the Ru-0 planes 
and U^(ij) = U^(ij) = U± for i and j being near- 
est on neighbouring planes explained most of the avail- 
able experimental data. Namely, the corresponding solu- 
tion featured a gap function on the 7-sheet of the form 
A cc {k) ~ smk x + isin/c y and and a line of gap on the 
(3 sheet. For this scenario the requirement that there be 
only one transition at T c ~ 1.5 K fixed both U±_ and 
U\\ and hence all further results could be regarded as 
quantitative predictions of the model. Remarkably, the 
model gave a satisfactory account of the data for the spe- 
cific heat C(T), superfluid density n s (T) and the thermal 
conductivity n(T). 

We have also investigated the stability of the model to 
introduction of further interaction constants and disor- 
der. We found that the predictions of the model are 
robust to changes of new interactions, while disorder 
mainly affects the f-wave solution. Thus we can conclude 
that the experimental data support a simple model which 
describes, quantitatively, the p-wave pairing observed in 
Sr2Ru04 on the basis of two orbital specific coupling con- 
stants: Uu = 40 meV U± = 48 meV. The central physical 
feature of the model is that Uu corresponds to interac- 
tion between electrons in the Ruthenium planes while 
U± describes an inter-plane interaction of roughly equal 
strength. 

In view of the above results, we would like to empha- 
size two points. Firstly, we have proposed an alternative 
to the the 'intra-band proximity effect' model of Zhito- 
mirsky and Rice 25 for describing horizontal line nodes on 
the a, P sheets of the Fermi Surface in superconducting 
Sr2Ru04. Our bond model differs from theirs in the way 
the interlayer coupling is implemented. The extension of 
the model in the spirit of ZR has also been studied by 
allowing for 3-site interactions in the Hamiltonian. Even 
though the resulting 'bond proximity model' features sin- 
gle superconducting transition temperature the original 
model with fine tuned two interactions gives better fit to 
experimental T dependence of the specific heat. 
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